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Abstract. One of the main problems in modeling the ionised outflows in Active Galactic Nuclei is the unknown distance of 
the outflowing wind to the central source. Only if the density is known this distance can be determined through the ionisation 
parameter. Here we study density diagnostics based upon O V transitions. O V is known to have metastable levels that are 
density dependent. We study the population of those levels under photoionisation equilibrium conditions and determine for 
which parameter range they can have a significant population. We find that resonance line trapping plays an important role in 
reducing the critical densities above which the metastable population becomes important. We investigate the K-shell absorption 
lines from these metastable levels. Provided that there is a sufficient population of the metastable levels, the corresponding 
K-shell absorption lines are detectable and are well separated from the main absorption line originating from the ground state. 
We then present the Chandra LETGS spectrum of the Seyfert 1 galaxy Mrk 279 that may show for the first time the presence 
of these metastable level absorption lines. A firm identification is not yet possible due to both uncertainties in the observed 
wavelength of the strongest line as well as uncertainties in the predicted wavelength. If the line is indeed due to absorption from 
O V, then we deduce a distance to the central source of one light week to a few light months, depending upon the importance 
of additional heating processes. 

Key words. Galaxies: individual: Mrk 279 - Galaxies: Seyfert - quasars: absorption lines - quasars: emission lines — X-rays: 
galaxies 

1 . Introduction ample of this broad range in ionisation parameter is the XMM- 

Newton Reflection Grating Spectrometer (RGS) spectrum of 
Active Galactic Nuclei (AGN) have X-ray spectra that are dom- NGC 5548i which shows the presence of K . shell absorption 

inated by nonthermal emission from the immediate surround- Unes from ^ oxyg£n i(ms between Q m _ Q ym (Steenbmgge 

ings of a supermassive black hole. Seyfert galaxies constitute f , fennel 
the lower luminosity class of AGN. In a significant fraction 

of all Seyfert galaxies the spectral signatures of the so-called ° ne of the main problems in modeling this outflow is 
warm absorber are visible (e.g. Reynolds fT997l George et al. that its distance r from the nucleus is hard to measure. 
1998 1. This warm absorber is very likely a photoionised out- Photoionisation modeling only yields the product nr 2 , as well 
flow from the accretion disk, with strong UV (for example as ^ column densit y nd with d the Sickness of the absorbing 
Crens haw et al. ll999l and X-ray absorption lines (Kaastra et al. la y er in the line of si S ht - Lacking a density, both the thick- 
120001 . X-ray observations with high sensitivity show the pres- ness _l and the distance r are essentially unknown. Gabel et 
ence of absorption lines with a broad range of ionisation pa- aL (2003) detected UV absorption lines from the metastable 
rameter £ = L/nr 2 , where L is the 1-1000 Ryd luminosity, n 2s 2 P 3p tti P let in C m aroymd 1 175 A in the hi S hest velocit y 
the density and r the distance to the central source. A good ex- outflow component in NGC 3783. These measurements indi- 
cate a relatively high density of order 10 15 m 3 correspond- 



Send offprint requests to: J.S. Kaastra ing to an upper limit of the distance to the nucleus of 0.3 pc. 

* Deceased 4 May 2004 However these density estimates are somewhat uncertain due 
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Galactic absorption in that source. Behar et al. 2003 pointed out 
that the 2s 2p triplet of C III can in fact be populated at low den- 
sities. Only one fine-structure level in the triplet requires high 



densities (see, e.g., Fig. 1 in Bhatia & Kastner 1993). Indeed 
updated calculations by Gabel et al. (2004a) show a lower den- 
sity of 3 x 10 10 m -3 and a distance of 25 pc. 

Other constraints on density hence distance may be de- 
rived from reverberation studies of the warm absorber in re- 
sponse to continuum variations, by measuring the recombina- 
tion time scale. Long grating observations of NGC 3783 with 
XMM-Newton (Behar et al. 2003 ) and with Chandra (Netzer 



et al. 2003) have allowed for attempts to study the response of 
the absorber to the ionising continuum through reverberation. 
However, since no response was detected these could provide 
only lower limits (of the order of 0.5 pc) to the distance of the 
absorber from the central source. Reeves et al. 2004 did claim 
to see variability in the Fe-K absorption of that source, albeit 
with CCD spectra. They deduce an upper limit to the distance 
of the absorber to the central source of 0.02 pc. 

In stellar coronae, with collisional ionisation equilibrium, 
several X-ray emission lines are density-dependent, for ex- 
ample the forbidden and intercombination lines of the O VII 
triplet. As these lines have low oscillator strengths, it is im- 
possible to observe them in absorption. In some isoelectronic 
sequences, close above the ground state of the ion there ex- 
ist metastable levels that can have a significant population. A 
well known example is the Be-sequence, where the C III line 
at 1909 A is often used as a density diagnostic in plasmas. 
Another ion in this iso-electronic sequence is O V, where the 
metastable ls 2 2s2p 3 P level is only 10 eV above the ground 
state ls 2 2s 2 1 So- The J = term of this triplet cannot decay 
radiatively to the ground state, leading to a significant popula- 
tion of this level at all densities. K-shell absorption from the 
ground state of this ion (at 22.3 A) has been observed in AGN 
(NGC 5548, Steenbrugge et al. (2003, 2004 1; NGC 4051, (Ogle 
et al. l2004> : and recently Mrk 279 (Costantini et al. l2004> . In 
this last source the spectrum indicates the possible presence 
of absorption from the metastable level of O V. The critical 
density for the population of this metastable level is around 
10 16 m -3 , a relevant density for AGN outflows. The presence 
or absence of K-shell absorption lines from the metastable level 
of O V then is an important density diagnostic, which can serve 
to constrain distance and geometry of the outflow. 

In principle transitions from the metastable level of O V 
can also be studied using the UV absorption lines from the 
2s 2p triplet to 2p 2 triplet near 760 A. Pettini & Boksenberg 
( 1986 1 identified these lines in IUE spectra of the BAL quasar 
PG 0946+301 . However, in HST data of the same object no ev- 
idence for these lines was found (Arav et al. ( 1999 1. But in yet 
another BAL quasar (QSO B0226-1024) Korista et al. ([1992) 
found evidence for the presence of these lines, implying den- 
sities of the order of 10 17 m -3 . However, due to Galactic ab- 
sorption such an analysis can only be done for bright, highly 
redshifted quasars. For nearby Seyfert galaxies K-shell X-ray 
absorption lines are the ideal tool to study densities. 

In this paper we study the expected population of the 
metastable level of O V under photoionised conditions, as well 
as the wavelengths and oscillator strengths of the correspond- 



ing K-shell absorption lines. We then present a possible de- 
tection of absorption from O V* levels in the Chandra LETGS 
spectrum of Mrk 279. 

2. Level populations 

2. 1. Transition rates 



Table 1. Energy levels of O V. Energies are taken from Wiese 
et al. ( 1996). Levels 2-4 are the metastable levels discussed in 
this paper. 



level 


Configuration 


J 


Energy (eV) 


1 


2s 2 *S 








2 


2s 2p 3 P 





10.16 


3 




1 


10.18 


4 




2 


10.21 


5 


2s2p!p 


1 


19.69 


6 


2p 23 P 





26.47 


7 




1 


26.49 


8 




2 


26.52 


9 


2p 2 X D 


2 


28.73 


10 


2p 21 S 





35.70 



We made a model for the population of all 10 levels with 
principal quantum number n = 2. These levels are listed in 
Tabled The energy-level diagram is shown in Fig.^ We omit 
the n > 2 levels from our calculation, as the lowest of these 
levels has an energy of 67.82 eV. This is too high as compared 
to the typical temperatures of a few eV for which O V is formed 
under photoionised conditions. 

The population levels n^ were solved using the equation 



ij ^j i 



(1) 



where we sum over j and Rij is the total transition rate matrix 
from level j to level i, and bi represents source/loss terms like 
ionisation from O IV to O V or recombination from O VI to O V. 
Other loss terms like ionisation or recombination of O V itself 
can be accomodated in R^, as these processes are proportional 
to the O V density. 

Radiative transition rates between the levels of Table[2were 
taken from the compilation of Wiese et al. ( 1996 1. 



We took the excitation rates from Safranova et al. (1995 1 
for the following transitions: 2s 2 -2s2p and 2s2p-2p 2 . The 
remaining rates were taken from Kato et al. (1990 1. Kato et 
al. give only the total rates for the following transitions in- 
volving a singlet and a triplet: 2s 2 -2p 2 3 P, 2s 2p 3 P-2s 2p 1 P, 
2p 2 3 P-2p 2 1 D and 2p 2 3 P-2p 2 1 S. Lacking more information, 
we therefore subdivided these rates according to the statisti- 
cal weights of the triplets involved. This is also found to be 
valid for the distorted wave calculations of Zhang & Sampson 



(1992 1, which explicitly give multiplet-resolved rates. There is 
a caveat, however. As noted by Safranova et al. ( 1995 1 in the 
forbidden transitions (J 7 — J = — 0,0 — 2,2 — 0), only the ex- 
change contribution to the collision strength is important, and 
the direct part is zero. The 3 P - 1 D 2 , 3 P - 1 S and 3 P 2 - 1 S 
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2p 2 2s2p2p 2 

1 S 1 P 1 D 
o Q r 1 u 2 



2s2p 



2p' 



3p 3p 3p 3p 3p 3p 
r r 1 r 2 r r 1 r 2 




Fig. 1. Energy level diagram of O V. Only the n = 2 levels 
are shown. The energy differences within the 3 P triplets have 
been exaggerated for clarity. Solid lines: transitions with tran- 
sition probabilities larger than Is -1 ; the line thickness is pro- 
portional to the logarithm of the transition probability. Dashed 
lines: transitions with transition probabilities less than 1 s . 



transitions of the above mentioned configurations are forbidden 
and would suffer from this effect. As however in most astro- 
physical situations the occupation of the 2p 2 levels is low, and 
these levels have strong radiative transitions to lower levels, the 
errors made by our approximation are not too serious. 

We added the proton excitation rates within the 2s 2p and 
2 p 2 3p tr ipi ets f rom Doyle et al. J1980I) . 



We also included radiative recombination (O VI to O V and 
O V to O IV) and collisional ionisation rates (O IV to O V and 
O V to O VI). We took the total radiative recombination rates 
from Arnaud & Rothenflug ( 1985 1 and assume that the rates 
towards the n — 2 levels, either direct or indirect, are propor- 
tional to the statistical weight of these levels. The collisional 
ionisation rates were also taken from Arnaud & Rothenflug 
(1985 1. For ionisation of Oiv, we assume that all OlV ions 
are in the 2s 2 2p 2 P ground state, and that 2s-shell ionisation 
of O IV results into a 2s 2p 3 P state of O V with the population 
of the sublevels distributed according to their statistical weight. 
Note however that for the parameter range that we consider the 



collisional ionisation and radiative recombination rates are rel- 
atively unimportant. 

At low temperatures charge transfer reactions with neutral 
hydrogen become important. We included the charge transfer 
recombination rates for recombination of O VI and O V from 
Rakovic et al. J2001E . Similar to the radiative recombination, 
we assume that the rates towards the n = 2 levels, either direct 
or indirect, are proportional to the statistical weight of these 
levels. 

2.2. Collisional ionisation equilibrium 

Using the rates determined in the previous section, we de- 
termine equilibrium population levels for a range of densi- 
ties and temperatures, first under collisional ionisation equi- 
librium (CIE) conditions. We restrict the temperature range to 
kT < 20 eV; higher temperatures are not expected under pho- 
toionised conditions, and moreover they would result in a sig- 
nificant amount of excitation to the n = 3 levels, with the cor- 
responding downwards radiative decays. 

For this temperature range, the effect of collisional ioni- 
sation is negligible, as the electrons have insufficient energy 
to overcome the 1 13.896 eV ionisation potential. Radiative re- 
combination affects the population of the 2s 2p 3 P triplet by 
no more than 3 %. Proton excitation modifies the populations 
only significantly (more than 3 %) for temperatures higher than 

10 eV, but even for kT — 30 eV the effects are smaller than 

11 % for all densities. 

Charge transfer recombination due to collisions with H I 
atoms appears to be very important at temperatures below 2 eV. 
This is mainly due to the large recombination cross section at 
low temperatures and the high neutral fraction at these low tem- 
peratures. Since under CIE conditions the Ov fraction is less 
than 0.001 for temperatures below 9 eV, inclusion of charge 
transfer recombination is of little practical use under those cir- 
cumstances, however. 

Fig-Elshows the relative occupation P of the 2s 2p 3 P triplet 
with respect to the ground state. The enhancement of this ratio 
at kT — 0.5 eV and kT = 1 eV is due to the charge transfer 
recombinations. The sharp increases in P as a function of hy- 
drogen density nn around densities of a) 10 12 , b) 10 17 and c) 
10 22 m -3 are due to the following. At a), the collisional cou- 
pling of the 2s 2p 3 P2 level with the 2s 2p 3 Pi level (either di- 
rect or via excitation to the 2p 2 3 Pi level followed by radiative 
decay) becomes stronger than the radiative decay by a forbid- 
den transition to the ground state. At b), the collisional coupling 
of the 2s 2p 3 Pi level with the other triplet levels is stronger 
than its radiative decay to the ground state. Finally, at c) the 
ground state depopulates due to strong collisional coupling of 
the 2s 2p *P level with the 2p 2 l T>, l S levels. This causes the ra- 
tio of the triplet with respect to the ground state to increase (the 
ratio with respect to the sum of all O V levels does not change). 

2.3. Photoionisation equilibrium 

For photoionised plasmas we use the same approach as before, 
only the transition matrix Rij and the source term bi contains 
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log n H (m ) 
16 




og n H (cm ) 

Fig. 2. Population of the 2s 2p 3 P triplet with respect to the 
ground state in O V as a function of density and kT (in eV) 
as labeled. All three sublevels ( J = 0, J = 1 and J = 2) have 
been added. 



more terms. We now include both photo-excitation and stim- 
ulated emission between the n — 2 levels. As we will use an 
ionising spectrum with a strong UV component (see below), 
photo-excitations towards higher levels than n = 2 are rela- 
tively unimportant. We include photoionisation of O IV from 
both the 2p and 2s shells. For our ionising spectrum, the pho- 
toionisation rate of O IV from the Is shell is only 2.5 % of the 
2p rate, and as Is ionisation of O IV is followed in most cases 
by an Auger transition, producing O VI, we neglect this process 
here. We also neglect the similar process of K-shell ionisation 
of O III followed by an Auger transition, as we assume that the 
concentration of O III relative to O V may be small and again 
K-shell photoionisation is rare relative to L-shell photoionisa- 
tion for our ionising spectrum. However, we included the 5 % 
contribution of Is photoionisation of Ov to the total n = 2 
photoionisation rate, as this represents a loss term for which 
we are not interested in the end product. 

For the photoionising spectrum we adopt a model based 
upon the observed Chandra and UV continuum of Mrk 279 
(Costantini et al. 2004). This model is similar to the model em- 
ployed for NGC 5548 by Kaastra et al. d2002l . The spectrum 
contains a strong UV bump (Fig. |3j. We normalized it to a 1- 
1000 Ryd luminosity L of 8.31 x 10 37 W, in order to reproduce 
the observed UV and X-ray continuum of Mrk 279. Note that 
for this spectral shape the relation between the photoionisation 
parameters £ and U is given by U — Q/4irr 2 nc — 0.036£, 
where Q is the number of hydrogen ionising photons and £ is 
expressed in units of 10~ 9 Wm. The concentrations of OlV- 
O VI as well as the temperature T(£) as a function of the ionisa- 
tion parameter £ were obtained from a set of runs with Cloudy 
Version 95.06 (Ferland 2002) for a thin photoionised slab at 
low density and with solar abundances. In that approximation, 
T(£) and the ion concentrations only depend upon the shape 
of the ionizing spectrum. Fig. |4] shows the Ov concentration 
and T(£) relation. We then determined the level populations 




100 1000 
Wavelength (k) 

Fig. 3. Spectral energy distribution (thick line) used for our 
photoionisation calculations. The thin line shows the same 
spectrum absorbed by a Galactic foreground of cold gas with 
a hydrogen column density of 1.64 x 10 24 m~ 2 and a v = 
lOkms- 1 . 




Fig. 4. Temperature (upper panel) and Ov fraction (lower 
panel) as a function of ionisation parameter £. The maximum 
Ov concentration of 40 % occurs for log£ = —0.05, at a 
temperature kT = 1.9 eV (to be compared with 47 % at 
kT = 21.0 eV under CIE conditions). 



as a function of «h and T; inverting T(£) yields £ and using 
£, L and nn we then know the distance of the absorber to the 
central source and hence the relevant fluxes for photoionisation 
and photo-excitation. 

We show our results in Fig. [5] At low temperatures (kT < 
1 eV) the neutral hydrogen fraction is now much smaller than 
in a collisional equilibrium plasma, resulting in less important 
effects of charge transfer recombination. For temperatures kT 
between 2-4 eV and densities below 10 20 m~ 3 the results for 
photoionisation equilibrium (PIE) and CIE are not significantly 
different. At higher densities photo-excitation of in particular 
the 2s 2 -2s2p *Pi resonance line (629.73 A) is more impor- 
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Fig. 5. Population of the 2s 2p 3 P triplet with respect to the 
ground state in O V as a function of density and kT (in e V) as 
labeled for a photoionised plasma. All three sublevels (J = 0, 
J = 1 and J = 2) have been added. 

tant than collisional excitation from the ground state to the 
2s 2p Pi level, causing enhanced higher level populations. 

2.4. Resonance line trapping 

In Sect. 12.31 our treatment of PIE assumed that O V is located 
in a photoionised slab with negligible column density. However 
modeling of the Mrk 279 spectrum shows that O V has a col- 
umn density of order 10 20 m~ 2 . For a velocity broadening o v 
of 50 km s _1 this implies that the 2s 2 -2s2sp 1 P 1 resonance line 
has an optical depth of 69 at the line center. As for a broad range 
of parameter space radiative decay through this emission line 
is the dominant decay proces of the upper level j, line pho- 
tons effectively undergo resonance line scattering. The direct 
escape probability P csc per emitted photon is therefore very 
small, and this leads effectively to a smaller radiative decay 
rate A',^ = P csc Aij, with A%j the true transition probability. 
We follow here the approach of Hollenbach & McKee ( 1979) 



and Kallman & McCray (1982i in order to approximate P ci 
for a slab: 



P™ 



(Aj/Atot) 



1. 



ToV / 27rln(2.13 + r 2 ) 



(2) 



where r is the optical depth at line center, and A tot is the to- 
tal transition probability from level j to any higher or lower 
level, including radiative and collisional processes. Resonance 
line trapping appears to be relevant for the transition between 
levels 1-5 (629.73 A) mentioned above, as well as for the 3-6 
(761.13 A), 5-9 (1371.30 A) and 5-10 (774.52 A) transitions 
whenever the 2s 2p configuration has a significant population. 
In Fig.[6]we show the results for our calculations including 
resonance line trapping. It is seen that the critical density at 
which the 2s 2p triplet level gets a significant population can 
shift downwards by several orders of magnitude for sufficiently 
large column densities. 
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Fig. 6. Population of the 2s 2p 3 P triplet with respect to the 
ground state in O V as a function of density and kT (in eV) in 
a photoionised plasma. Solid curves are for kT = 2 eV, dash- 
dotted curves for kT — 4 eV. Labels indicate logiV with N 
the column density of O V in m~ 2 . All three sublevels ( J = 0, 
J = 1 and J = 2) have been added. 



3. X-ray absorption lines 

We have calculated wavelengths and oscillator strengths for all 
transitions between levels with a K-shell vacancy (ls2s 2 2p, 
Is2s2p 2 and ls2p 3 ) and the Is 2 2s 2 , ls 2 2s2p, ls 2 2p 2 and 
ls 2 2snp levels. We used the code of Cowan (1981) for this 
purpose. 

Table [2] lists the calculated wavelenghts and oscillator 
strengths / of all lines with / > 0.001 and principal quantum 
number n < 10. We determined only total oscillator strengths 
and average wavelengths of the multiplets, as the current X-ray 
instrumentation has insufficient resolution to resolve the indi- 
vidual lines of the multiplets. 

An important question is the accuracy of the wavelengths. 
In general these wavelengths have never been measured in the 
laboratory, with the exception of the resonance transition A2. 
This wavelength is 22.374±0.003 A (Schmidt et al. l2004l 

The lower levels (Is 2 2£2£') of these transitions have ac- 
curately known energy levels Ei (Wiese et al. 1996 1. For the 
upper levels ( 1 s 21 21' n£") accurate energy levels Ej can be 
obtained indirectly from Auger electron spectra through the re- 
lation Ej = Ea + I, with I the ionisation potential of O V, 
for which we take 113.896 eV (Wiese et al.Q]J96), and E A the 
energy of the Auger electron. 

Ea has been measured in the laboratory or calculated the- 
oretically; we list the relevant values in Table [3] A compari- 
son of the calculations of the group of K.T Chung (Lin et al. 
2001 Shiu et al. 2001 Lin et al. 2002 1 with the measurements 
of Bruch et al. ( 119791 fl987i shows excellent agreement; the 
average difference between observed and calculated values is 
—0.05 ±0.1 6 eV; we therefore adopt a nominal uncertainty of 
0.20 eV for these calculated energies, equal to the measurement 
errors in the best measurements. A comparison of the energies 
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Table 2. Strongest X-ray lines of O V. The values listed apply 
for the full multiplets (not splitted into sublevels). 



label 



Initial 
conf 



state 
term 



Final 



conf 



state 
term 



/ 



A 

(A) 



A2 


2s 2 


*S 


s2s 2 2p 


X P 


0.649 


22.381 


A3 


2s 2 


: S 1 


s2s 2 3p 


l V 


0.108 


19.871 


A4 


2s- 


: S ] 


s 2s 2 4p 


*P 


0.041 


19.258 


A5 


2s- 


: S 1 


s2s 2 5p 


*P 


0.020 


19.002 


A6 


2s 2 


: S 1 


s2s 2 6p 


X P 


0.012 


18.869 


A7 


2s 2 


: S ] 


s2s 2 7p 


J P 


0.007 


18.791 


A8 


2s 2 


: S 1 


s2s 2 8p 


X P 


0.005 


18.742 


A9 


2s 2 


: S ] 


s2s 2 9p 


*P 


0.004 


18.708 


Bl 


2s 2p 


3p j 


s 2s 2p 2 


( 3 S) 3 P 


0.328 


22.466 


B2 


2s 2p 


3p ] 


s 2s 2p 2 


( 3 S) 3 D 


0.178 


22.431 


B3 


2s 2p 


3p j 


s 2s 2p 2 


( 3 S) 3 S 


0.039 


22.239 


B4 


2s 2p 


3p j 


s 2s 2p 2 


( X S) 3 P 


0.014 


22.132 


CI 


2s 2p 


: P ] 


s 2s 2p 2 


(^D 


0.182 


22.558 


C2 


2s 2p 


*P 1 


s2s2p 2 


(?$fv 


0.339 


22.389 


C3 


2s 2p 


X P 


s2s2p 2 


^S^S 


0.040 


22.363 


Dl 


2p 2 


3p 


s 2s 2 2p 


3p 


0.005 


23.643 


D2 


2 P 2 


3p 


s2p 3 


! D 


0.173 


22.511 


D3 


2 P 2 


3p 


s2p 3 


3 S 


0.151 


22.474 


D4 


2 P 2 


3p 


s2p 3 


3p 


0.105 


22.339 


El 


2r 


: D 


s 2s 2 2p 


l P 


0.006 


23.598 


E2 


2p 2 


: D 


s2p 3 


l D 


0.328 


22.477 


E3 


2r 


: D 


s2p 3 


3p 


0.002 


22.457 


E4 


2 P 2 


: D 


s2p 3 


l P 


0.112 


22.306 


Fl 


2 P 2 


: s 


s 2s 2 2p 


l P 


0.002 


23.855 


F2 


2p 2 


: s 


s2p 3 


X P 


0.438 


22.536 



calculated by Chen (1985 1 with the above mentioned calcula- 
tions by the group of K.T Chung shows an average difference 
of 1.08±0.28 eV, with a standard deviation of 0.60 eV for the 
9 levels in common. Therefore whenever we use the energies 
calculated by Chen, we subtract 1.08 eV and assign a nominal 
accuracy of 0.60 eV. 

The adopted energies Ej of Table [3] were combined with 
the lower level energies Ei of Table ^ m order to calculate 
the predicted wavelengths of the transitions from Table |2] We 
list those wavelengths in Table |4] Champeaux et al. (2003) 
recently deduced the ionisation potential / for Ov to be 
113.66±0.13 eV, based upon the quantum defect estimated 
from resonances in the measured photoionisation cross section. 
If this value is used instead of the value of 113.896 eV given 



by Wiese et al. ( 1996 1, most wavelengths listed in the column 
"Auger" of TableHincrease by 0.009 A. 

A comparison of the energies of the ls 2 2£2£' levels cal- 
culated by us using the Cowan (1981 1 code with the values 
of Table [2 (Wiese et al. 1996l shows that our calculated ener- 
gies are smaller by 0.71 (levels 2-4), 0.91 (level 5), 0.83 (lev- 
els 6-8), 0.16 (level 9) and 1.46 eV (level 10), respectively. 
Correcting all wavelengths where these levels are involved by 
these differences and also adjusting all K-shell vacancy states 
by the same number so as to get the laboratory wavelength of 
the 22.374 A line correct, yields the column labeled with "corr." 
in Table|4] The difference in wavelength between both columns 
can be as large as 0.053 A, comparable to the energy resolution 
of Chandra's LETGS. Finally, we list in Table |4] also wave- 



Table 3. Core-excited energy levels in Ov. Columns (1) and 
(2): configuration and term; (3) and (4): Calculated energy 
of the Auger electron E& for a transition to the Is 2 2s 2 S 
(ground) state of O VI; where only transition energies to the 
Is 2 2p 2 P state are given in the original reference these have 
been corrected for the energy difference between the 2 S and 2 P 
state; (5): measured Ea', (6): adopted Ea', (7): adopted energy 
Ej = E A + I of the level. 



conf. 



term 



E A (th) 
(eV) 



E A (th) 
(eV) 



E A (exp) 
(eV) 



E A 
(eV) 



(eV) 



ls2s 2 2p X P 


441.16" 


440.55° 


440.5±0.2 / 


440.55 554.45±0.20 


3p 


436.63 a 


435.96 6 


435.9±0.2 / 


435.96 549.86±0.20 


ls2s 2 3p X P 




507.54° 




507.54 621.44±0.20 


Is2s2p 2 X S 


462.10 a 


460.83 c 


460±2 e 


460.83 574.73±0.20 


*P 


461. 98 a 


460. 12 c 


460.5±T 


460.12 574.02±0.20 


*D 


456.74 a 


455.35 c 


456.3±0.8 e 


455.35 569.25±0.20 


3 S 


455.06° 


454. 17 b 


454±l e 


454.17 568.07±0.20 


( 3 S) 3 P 448.76° 


448. 70 6 


448.6±0.2 / 


448.70 562.60±0.20 


^SfP 458.58" 


456.9 e 




457.50 571.40±0.60 


3 D 


449.74° 


448.68 b 


448.5±l e 


448.68 562.58±0.20 


5p 


438.33° 




438.55±0.2 J 


438.55 552.45±0.20 


ls2p 3 X P 


473.72° 


471.80° 




471.80 585.70±0.20 


*D 


468.05° 






466.97 580.87±0.60 


3 S 


464.91° 






463.83 577.73±0.60 


3p 


469.40° 






468.32 582.22±0.60 


3 D 


463.75° 






462.67 576.57±0.60 


5 S 


456.27° 




456.4±0.2 / 


456.40 570.30±0.20 


a ChenfT985l 








b Linetal. 


20011 








c Shiu et a 


■I200T1 








d Lin et al.EQOa 








c Bruch et al.lT979l 








1 Bruch et al.lT987l 









lengths from the paper by Chen &Crasemann( 1987 1 as well as 
wavelengths presently calculated with the Hebrew University 
Lawrence Livermore Atomic Code (HULL AC). 

In general, the theoretical calculations for the wavelengths 
differ quite a lot. We illustrate this with the resonance line at 
22.374 A. Calculated values range from 22.334 A (HULLAC) 
to 22.415 A (Chen & Crasemann); for comparison, Pradhan et 
al. (I20"0"3l give a value of 22.35 A. 

There is a good agreement between our calculations for os- 
cillator strengths and transition probabilities and those of Chen 
& Crasemann ( 1987 1, see Table|5] The rms difference between 
both sets of oscillator strength is 0.003. Also, the total radiative 
transition probabilities from core-excited levels agree within 
2 %. It is necessary to take also the Auger transitions into ac- 
count for estimating the lifetime of the core-excited levels. We 
use the Auger rates determined by Chen & Crasemann ( 1987 1. 
Their rates are in good agreement with the calculations of Lin et 
al. J200 1 1 120021 for the 6 levels in common. These Auger rates 
are an order of magnitude larger than the radiative transition 
rates. Thus neglecting Auger rates will affect the estimated line 
equivalent width for high column densities. There are no Auger 
rates available for the A4-A9 transitions. Lin et al. (2002 1 cal- 
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Table 4. Wavelength comparison of the X-ray lines of O V. 
Wavelengths are given in A. The column "Cowan" gives A from 
Table |2j "corr." are these calculated wavelengths corrected as 
described in the text; "Auger" are the wavelengths determined 
from Auger energies as listed in Table [5] "CC" are the cal- 



culations by Chen & Crasemann (1987); "HULLAC" are the 
wavelengths calculated using the HULLAC code. 



line Cowan corr. 



A2 
A3 
A4 
A5 
A6 
A7 
A8 
A9 
Bl 
B2 
B3 
B4 
CI 
C2 
C3 
Dl 
D2 
D3 
D4 
El 
E2 
E3 
E4 
Fl 
F2 



Auger 



CC 



22.381 
19.871 
19.258 
19.002 
18.869 
18.791 
18.742 
18.708 
22.466 
22.431 
22.239 
22.132 
22.558 
22.389 
22.363 
23.643 
22.511 
22.474 
22.339 
23.598 
22.477 
22.457 
22.306 
23.855 
22.536 



22.374 
19.866 
19.253 
18.997 
18.864 
18.786 
18.737 
18.703 
22.488 
22.453 
22.260 
22.153 
22.589 
22.419 
22.393 
23.670 
22.538 
22.501 
22.365 
23.597 
22.476 
22.456 
22.305 
23.908 
22.589 



22.362± 0.008 
19.951± 0.008 



22.415 



22.444± 
22.445± 
22.224± 
22.092± 
22.561± 
22.367± 
22.338± 
23.690± 
22.540± 
22.492± 
22.311± 
23.584± 
22.455± 
22.400± 
22.260± 
23.901± 
22.543± 



0.008 
0.008 
0.008 
0.024 
0.008 
0.008 
0.008 
0.009 
0.024 
0.024 
0.024 
0.009 
0.024 
0.024 
0.008 
0.009 
0.008 



22.535 
22.494 
22.279 
22.139 
22.668 
22.453 
22.448 
23.790 
22.615 
22.566 
22.383 
23.747 
22.584 
22.529 
22.353 
24.112 
22.677 



HULLAC 



22.372 
19.921 
19.316 
19.060 
18.949 
18.849 
18.800 
18.766 
22.425 
22.415 
22.197 
22.091 
22.522 
22.365 
22.349 
23.572 
22.478 
22.404 
22.331 
23.520 
22.430 
22.479 
22.283 
23.765 
22.503 



Table 5. Oscillator strengths /, total transition probabilities 
A ra d, Auger decay rates Aa and total transition rates A to t of 
O V. Transition probabilities are in units of 10 12 s _1 . 



id 



r 



^ra 



A A a " 



a Chen & Crasemann fT987l 

b Present work 

c for A3-A9 see text 



A tl 



A2 


0.646 


0.649 


3.05 


3.00 


87.9 


90.9 


A3 




0.108 




0.63 


129 


129.6 


A4 




0.041 




0.26 


129 


129.3 


A5 




0.020 




0.13 


129 


129.1 


A6 




0.012 




0.07 


129 


129.1 


A7 




0.007 




0.05 


129 


129.1 


A8 




0.005 




0.03 


129 


129.0 


A9 




0.004 




0.02 


129 


129.0 


Bl 


0.326 


0.328 


4.29 


4.34 


34.2 


38.5 


B2 


0.171 


0.178 


1.35 


1.41 


94.2 


95.6 


B3 


0.039 


0.039 


1.56 


1.56 


46.9 


48.5 


B4 


0.014 


0.014 


0.19 


0.19 


100 


100.2 


CI 


0.181 


0.182 


1.41 


1.43 


213 


214.4 


C2 


0.333 


0.339 


4.40 


4.51 


22.9 


27.4 


C3 


0.041 


0.040 


1.63 


1.60 


154 


155.6 


Dl 


0.007 


0.005 


0.08 


0.06 


129 


129.1 


D2 


0.169 


0.173 


1.32 


1.37 


131 


132.4 


D3 


0.153 


0.151 


6.00 


5.99 


0.0015 


6.0 


D4 


0.102 


0.105 


1.36 


1.42 


79.9 


81.3 


El 


0.009 


0.006 


3.05 


3.00 


87.9 


90.9 


E2 


0.330 


0.328 


4.31 


4.35 


130 


134.4 


E3 


0.000 


0.002 


1.36 


1.42 


79.9 


81.3 


E4 


0.111 


0.112 


4.36 


4.41 


76.9 


81.3 


Fl 


0.003 


0.002 


3.05 


3.00 


87.9 


90.9 


F2 


0.439 


0.438 


4.36 


4.41 


76.9 


81.3 



Using A t 



A ra d + Aa with j4 ra d from the present work 



culate an Auger rate of 1.29 x 10 s _1 for the upper level 
of transition A3. In their calculations for Be, they find that for 
higher principal quantum number n the Auger rates do not de- 
pend strongly upon n. Therefore we adopt a constant value for 
A a for transitions A3-A9. 



4. Comparison with the observed spectrum of 
Mrk 279 

4. 1. Chandra LETGS observations 

Mrk 279 is a moderately distant (z — 0.0306) Seyfert 1 galaxy 
that we observed for 340 ks with the Chandra LETGS. This in- 
strument has a spectral resolution of ^0.05 A FWHM. The full 
data analysis of the LETGS spectrum of Mrk 279 will be pre- 
sented elsewhere (Costantini et al. 2004 1. Here we only sum- 
marize the results relevant for the study of the O V absorp- 
tion lines. The LETGS spectrum was fit by a simple power 
law plus modified blackbody with Galactic absorption. This 
smooth continuum spectrum produces an excellent fit to the 
data except for the few regions where Mrk 279 has absorp- 
tion lines. The column densities in Mrk 279 are several times 
smaller than in other well-studied AGN, and therefore most ab- 
sorption lines and edges are very weak except for the lines from 



the most abundant metal, which is oxygen. These lines are pre- 
dominantly found in the 18-24 A band. 

Our continuum fit is excellent everywhere except for the 
region near the K-shell transitions of oxygen (Fig. [7}. Here a 
broad emission feature with a peak amplitude of about 30 % of 
the continuum is visible. Similar broad X-ray emission lines 
have been seen before in other AGN (Kaastra et al. 2002 
Steenbrugge et al. 2004i, and may have asymmetric line pro- 
files. This broad line is due to the O VII triplet from the broad 
line region. The presence of this broad feature makes the es- 
timate of the underlying continuum for the narrow absorption 
lines less certain. We estimated the broad line profile by draw- 
ing a spline through the spectral regions that are known to be 
line-free. 

From a detailed study of the Ovi UV and X-ray lines 
(Gabel et al. 2004b, Costantini et al. I2004J we find that the 
dominant absorber has an outflow velocity of 300 km/s and 
a line profile that can be approximated to first order by a 
Gaussian with ct v = 50 kms -1 . 
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21.5 



22 22.5 

Mrk 279 restfrome wavelength (A) 



Fig. 7. Chandra LETGS spectrum of Mrk 279 near the oxy- 
gen absorption lines. The observed spectrum has been divided 
by the best fit simple continuum model described in the text 
(dashed line). The additional broad emission feature is indi- 
cated by the curved line. The wavelength scale is for the rest- 
frame of Mrk 279 (at z = 0.0306). In the line identifica- 
tions, we have taken into account a small outflow (blueshift) 
of 300 km/s with respect to this restframe; this outflow has 
been measured using simultaneous high-resolution UV spec- 
tra. A few Galactic absorption lines are labeled with "Gal". 



4.2. Observed O v lines 

We clearly see the absorption at 22.37 A, which we identify 
as the Is 2 2s 2 - ls2s 2 2p 1 P 1 transition of Ov (Fig. Q. Its 
measured wavelength in the rest frame of Mrk 279, corrected 
for the 300 km/s outflow velocity, is 22.377±0.009 A, in ex- 
cellent agreement with the laboratory wavelength of 22.374 A 
mentioned before (Schmidt et al. 2004 The quoted error bar 
is only the statistical error; we estimate the systematic uncer- 
tainty (due to calibration uncertainties in the HRC-S detector) 
to be 0.008 A. 



Table 6. Observed absorption lines in Mrk 279, due to or re- 
lated to O V. The observed wavelength A is given with two er- 
ror bars: the first is the statistical uncertainty, the second the 
systematic due to the wavelength calibration uncertainty of the 
LETGS. For lines with wavelength fixed, the observed equiv- 
alent width (EW) is merely the formal upper limit. The line 
around 22.50 A is either fitted by a single line (22.540 A) or 
alternatively by two lines (22.501 and 22.548 A). 



A (A) 


EW(rr 


A) 


identification 


22.377±0.009±0.008 


18±6 




A2 


19.957 (fixed) 


2±8 




A3 


19.325 (fixed) 


6±6 




A4 


22.540±0.009±0.010 


20±5 




B1+B2?; see text 


22.501±0.016±0.010 


10±6 




B1+B2?; see text 


22.548±0.010±0.010 


18±6 




?; see text 


22.269±0.015±0.006 


10±6 




B3? ; see text 




20 
log N v (m~ 2 ) 

Fig. 8. Curve of growth for the most important K-shell absorp- 
tion lines of O V, for a Gaussian velocity broadening a v of 
SOkms- 1 . 



We find no clear evidence for transition A3. Upper limits 
to its equivalent width depend upon the adopted wavelength. In 
another Seyfert 1 galaxy with a strong A2 line we find evidence 
for a 2er detection of transition A3 at a laboratory wavelength of 
19.957±0.015 A, right within the wavelength range indicated 
in Table 0] We adopt this wavelength for transition A3 and use 
it to measure the limits to the line equivalent width (Table [6). 
We also do not detect line A4. We take here the wavelength as 
calculated with the HULLAC code to estimate the equivalent 
width. Using the velocity broadening er v of 50 km/s and the 
equivalent widths of lines A2, A3 and A4, we find a column 



density of 10 



20.16±0.42 „,-2 



m for O V in its ground state. The 



large uncertainty is mainly due to the strong saturation of the 
A2 line, which has an optical depth of 3.7 at its line center. 

Most absorption lines in the spectral band displayed in 
Fig-Elcan be easily identified by Galactic absorption lines or 
absorption lines due to the warm absorber. Only the feature 
around 22.50 A has no obvious identification with absorption 
from the ground state of any abundant ion, either in the warm 
absorber or the Galactic foreground. When fitted by a single 
line, its laboratory wavelength, corrected for the 300 km/s out- 
flow, is 22.540±0.009 A. A marginally better fit (A X 2 = -3.6 
for two additional degrees of freedom) is obtained for a double, 
blended line, with laboratory wavelengths of 22.501±0.016 A 
and 22. 548 ±0.010 A, respectively (see Table|6j. Again, quoted 
errors are statistical only; the systematic uncertainty is 0.010 A 
and is not correlated with the systematic uncertainty in the main 
22.374 A line. 

We suggest that in the case of a double blended line, the 
22.501 A component of the blend may be due to a combination 
of the Bl and B2 lines of O V, which have a wavelength in the 
22.47-22.53 A range. For the single line solution (22.540 A) 
this is somewhat less likely. We note that the added statisti- 
cal and systematic uncertainty of the measurement of the line 
centroid is in the 0.02-0.03 A range. The problem with a firm 
identification lies of course in both the uncertainty in the theo- 
retical wavelength as well as in the measured wavelength. 
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If our identification is true, we can deduce the column den- 
sity of O V in its 2s 2p 3 P triplet state. From the curve of growth 
(see Hig. |8j and an equivalent width of 10±6 mA, we find a 
column of io 19 - 82±0 - 44 m -2 . This combined with the column 
density of O V in its ground state yields a value for the popu- 
lation of the 2s 2p 3 P triplet state relative to the ground state 
of 0.5, with rms error bounds of a factor of 4 (i.e, the relative 
population is between 0.125 and 2). From the column density 
of the 2s 2p 3 P triplet state we predict that line B3 should have 
an equivalent width of 1 mA with an upper limit of 3 mA. The 
feature at 22.269 A (see Table|6jl could be marginally consistent 
with this. 

Table list the densities derived from the population ratio 
for a few temperatures. 

Table 7. Densities (in m~ 3 ) for 2s 2p triplet to ground ratio's 
of 0.125, 0.5 and 2, for a total Ov column density of 1.5 x 
m" 



kT(eV) 0.125 



0.5 



2.5 x 10 

10 17 
5 x 10 1 



TfT" 



3 x 10 
10 
4x 10 1 



TTT 



1N 



00 

1.5 x 10 19 
10 19 



From our photoionisation equilibrium calculations, it fol- 
lows that Ov has its maximum concentration for kT = 
2 eV (see Fig. |4j. This yields a lower limit to the density 
of 2.5 x 10 18 m~ 3 . kT = 2 eV corresponds to £ = 0.9 x 
10~ 9 Wm -2 m 3 , which then leads to an upper limit to the dis- 
tance r of 2 x 10 14 m or one light week. This would place the 
O V absorber close to the inner broad line region. Although not 
completely impossible, such a small radius is not very plausi- 
ble, because Balmer line time-delay measurements suggest a 
broad emission line region size of 6-17 light days for Mrk 279 
(Scott et al. 2004 and references therein), and UV spectra show 
that the broad lines are covered by the warm absorber. 

The temperature may be higher than 2 eV, however, if 
additional heating processes occur in the absorbing region. 
Cloudy (Ferland 2002 1 allows for the inclusion of additional 
heat sources. A temperature of 6 eV is reached for £ = 
0.7 x 10~ 9 Wm _2 m 3 , with an additional heating rate of 
1.0 x 10 _7 nf 4 Wm~ 3 , where nu the density in units of 
10 14 m~ 3 . In this case the additional heat source dominates 
the total heating rate. This is also the case for all temperatures 
above 3 eV, where the additional heating contributes more than 
80 % to the total heating rate. We find for kT = 6 eV a lower 
limit to the density of 5 x 10 16 m -3 , corresponding to an up- 
per limit to the distance of 1.5 x 10 15 m or 59 light days. This 
would place the absorber outside the broad line region. 

Is there a possible heating process that can produce the ad- 
ditional heat? For the lower limit to the density of 5 x 10 16 m~ 3 , 
the required heating rate Q is 0.025 Wm" 3 , One possibility 
would be classical heat conduction from a neighbouring hot re- 
gion. Using a hydrogen column density of 8 x 10 23 m~ 2 (cor- 
responding to a typical O V column of 1.5 x 10 20 m~ 2 ) and the 
lower limit of 5 x 10 16 m~ 3 for the density, we find an upper 
limit to the size d of the absorber of 1.6 x 10 7 m. For a heat 



conductivity k of 10" 10 T^ /2 WK" 1 m -1 with T K the temper- 
ature in K, and putting the heating rate Q due to conduction 
equal to nT/d 2 , we find that we need a temperature of 300 eV 
for the hot gas surrounding the absorbing medium. Hotter gas 
in photoionisation equilibrium coexisting with colder gas is in- 
deed known to be possible. As the geometry of the absorbing 
regions is still not fully understood, we leave this scenario as 
an interesting possibility. 

In our above analysis, we assumed that as far as O V is con- 
cerned the physical conditions in the absorber are such that O V 
has its maximum relative concentration. It is of course also pos- 
sible to have a hotter absorber with a smaller concentration of 
Ov. For instance, without invoking an additional heat source 
we find that kT = 4 eV if £ is 20 times larger than we used be- 
fore, and correspondingly the O v/O ratio is 100 times smaller. 
In that case most of the oxygen should be in the form of O VII. 
Our column densities for O VI and O VII (see Costantini et al. 



2004) are not inconsistent with such a scenario. In this case, 
we avoid the additional heating problem, but of course we still 
require a high density n > 10 17 m -3 in order to get the high 
metastable population fraction. In this case however the dis- 
tance to the central source is only 2 x 10 14 m or 8 light days. 

We do not expect blending due to lines from higher levels 
to be important for all temperatures below 4 eV (the line series 
labeled with C-F of Table|2ji. This is because for those temper- 
atures the population of these higher levels is small. 

Finally, if the metastable level has a significant popula- 
tion, also the UV absorption lines from this level should be 
detectable in principle. Unfortunately the strong, allowed UV 
absorption lines from the metastable level are in the unobserved 
part of the UV spectrum (near 760 A). If the density would be 
very high such that also the 2s 2p : Pi level gets a high popula- 
tion, then the strongest UV absorption line should be the transi- 
tion between levels 5-9 of Table[Jat 1371.30 A. Simultaneous 
HST observations of Mrk 279 (Gabel et al. 2004b I yield an 
upper limit to the column density for level 5 deduced from 
the non-detection of this line of 5 x 10 16 m~ 2 . Combined 
with the allowed range for the ground state column density 
(10 1982±044 m~ 2 as quoted before) this leads to an upper limit 
of 10~ 3 for the population of level 5 with respect to the ground. 
The corresponding upper limit to the density of 5 x 10 17 m~ 3 is 
consistent with the lower limits derived from the X-ray transi- 
tion from the metastable level. It should be noted, however, that 
the X-ray derived density constraints depend on several mod- 
eling and observational factors, as is evident from our present 
discussion. 

5. Conclusions 

We have investigated in this paper the population of the 
metastable level of O V under photoionised conditions. For suf- 
ficiently high densities and a temperature of just a few eV 
the metastable levels can have a significant population. The 
relevant critical densities also depend strongly upon the col- 
umn density, due to resonance line trapping. The population 
of higher levels than the 2s 2p 3 P triplet is always relatively 
small. If the metastable level has a significant population, then 
we should observe K-shell absorption lines from these levels. 
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The wavelengths of these lines are not very accurately known, 
however. This makes the identification of these lines in the X- 
ray spectrum of the Seyfert 1 galaxy Mrk 279 somewhat uncer- 
tain. If present, however, then the observed equivalent widths 
imply densities of order 10 17 m~ 3 or more, corresponding to a 
distance to the central source of the order of a light month. 
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